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Glossary 

• Black Box Model. This model assumes we can collect knowledge about an input / 
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thru queries without knowing how the answer to the query is computed. A synonym 
for black box is oracle. 

• Classical Computer. A computer which does not use the principles of quantum 
computing to carry out its computations. 

• Computational Complexity. In this article, complexity for brevity. The minimal 
cost of solving a problem by an algorithm. Some authors use the word complexity 
when cost would be preferable. An upper bound on the complexity is given by the cost 
of an algorithm. A lower bound is given by a theorem which states there cannot be an 
algorithm which does better. 

• Continuous Problem. A problem involving real or complex functions of real or 
complex variables. Examples of continuous problem are integrals, path integrals, and 
partial differential equations. 

• Cost of an Algorithm. The price of executing an algorithm. The cost depends on 
the model of computation. 

• Discrete Problem. A problem whose inputs are from a countable set. Examples of 
discrete problems are integer factorization, traveling salesman and satisfiability. 

• e- Approximation. Most real-world continuous problems can only be solved numeri- 
cally and therefore approximately, that is to within an error threshold e. The definition 
of ^-approximation depends on the setting. See worst-case setting, randomized setting, 
quantum setting. 

• Information-Based Complexity. The discipline that studies algorithms and com- 
plexity of continuous problems. 

• Model of Computation. The rules stating what is permitted in a computation and 
how much it costs. The model of computation is an abstraction of a physical computer. 
Examples of models are Turing machines, real number model, quantum circuit model. 

• Optimal Algorithm. An algorithm whose cost equals the complexity of the problem. 

• Promise. A statement of what is known about a problem a priori before any queries 
are made. An example in quantum computation is the promise that an unknown 1-bit 
function is constant or balanced. In information-based complexity a promise is also 
called global information. 

• Quantum Computing Speedup. The amount by which a quantum computer can 
solve a problem faster than a classical computer. To compute the speedup one must 
know the classical complexity and it's desirable to also know the quantum complexity. 
Grover proved quadratic speedup for search in an unstructured database. Its only con- 
jectured that Shor's algorithm provides exponential speedup for integer factorization 
since the classical complexity is unknown. 
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• Query. One obtains knowledge about a particular input thru queries. For example, 
if the problem is numerical approximation of f Q f(x) dx a query might be the evalu- 
ation of / at a point. In information-based complexity the same concept is called an 
information operation. 

• Quantum Setting. There are a number of quantum settings. An example is a 
guarantee of error at most e with probability greater than \. 

• Qubit Complexity. The minimal number of qubits to solve a problem. 

• Query Complexity. The minimal number of queries required to solve the problem. 

• Randomized Setting. In this setting the expected error with respect to the prob- 
ability measure generating the random variables is at most e. The computation is 
randomized. An important example of a randomized algorithm is the Monte Carlo 
method. 

• Worst-Case Setting. In this setting an error of at most e is guaranteed for all inputs 
satisfying the promise. The computation is deterministic. 

I Abstract 

Most continuous mathematical formulations arising in science and engineering can only be 
solved numerically and therefore approximately. We shall always assume that we're dealing 
with a numerical approximation to the solution. 

There are two major motivations for studying quantum algorithms and complexity for 
continuous problems. 

1. Are quantum computers more powerful than classical computers for important scientific 
problems? How much more powerful? 

2. Many important scientific and engineering problems have continuous formulations. 
These problems occur in fields such as physics, chemistry, engineering and finance. 
The continuous formulations include path integration, partial differential equations (in 
particular, the Schrodinger equation) and continuous optimization. 

To answer the first question we must know the classical computational complexity (for 
brevity, complexity) of the problem. There have been decades of research on the classical 
complexity of continuous problems in the field of information-based complexity. The reason 
we know the complexity of many continuous problems is that we can use adversary arguments 
to obtain their query complexity. This may be contrasted with the classical complexity of 
discrete problems where we have only conjectures such as P 7^ NP. Even the classical com- 
plexity of the factorization of large integers is unknown. Knowing the classical complexity of 
a continuous problem we obtain the quantum computation speedup if we know the quantum 
complexity. If we know an upper bound on the quantum complexity through the cost of a 
particular quantum algorithm then we can obtain a lower bound on the quantum speedup. 
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Regarding the second motivation, in this article we'll report on high-dimensional inte- 
gration, path integration, Feynman path integration, the smallest eigenvalue of a differential 
equation, approximation, partial differential equations, ordinary differential equations and 
gradient estimation. We'll also briefly report on the simulation of quantum systems on a 
quantum computer. 

II Introduction 

We provide a summary of the contents of the article. 

Section HTT1 Overview of quantum algorithms and complexity. 

We define basic concepts and notation including quantum algorithm, continuous problem, 
query complexity and qubit complexity. 

Section HV1 Integration 

High-dimensional integration, often in hundreds or thousands of variables, is one of the 
most commonly occurring continuous problems in science. In Section IIV.II we report on 
complexity results on a classical computer. For illustration we begin with a one-dimensional 
problem and give a simple example of an adversary argument. We then move on the d- 
dimensional case and indicate that in the worst case the complexity is exponential in d; the 
problem suffers the curse of dimensionality. The curse can be broken by the Monte Carlo 
algorithm. In Section HV.2I we report on the algorithms and complexity results on a quantum 
computer. Under certain assumptions on the promise the quantum query complexity enjoys 
exponential speedup over classical worst case query complexity. 

A number of the problems we'll discuss enjoy exponential quantum query speedup. This 
does not contradict Beals et al. [6] who prove that speedup can only be polynomial. The 
reason is that [6] deals with problems concerning total Boolean functions. 

For many classes of integrands there is quadratic speedup over the classical randomized 
query complexity. This is the same speedup as enjoyed by Grover's search algorithm of an 
unstructured database. To obtain the quantum query complexity one has to give matching 
upper and lower bounds. As usual the upper bound is given by an algorithm, the lower 
bound by a theorem. The upper bound is given by the amplitude amplification algorithm 
of Brassard et al. [12]. We outline a method for approximating high-dimensional integrals 
using this algorithm. The quantum query complexity lower bounds for integration are based 
on the lower bounds of Nayak and Wu [17] for computing the mean of a Boolean function. 

Section [V] Path Integration 

We define a path integral and provide an example due to Feynman. In Section IV. II we 
report on complexity results on a classical computer. If the promise is that the class of 
integrands has finite smoothness, then path integration is intractable in the worst case. If 
the promise is that the class of integrands consists of entire functions the query complexity 
is tractable even in the worst case. For smooth functions intractability is broken by Monte 
Carlo. In Section IV. 21 we report on the algorithm and complexity on a quantum computer. 
The quantum query complexity enjoys Grover-type speedup over classical randomized query 
complexity. We outline the quantum algorithm. 
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Section [VTl Feynman-Kac path integration. 

The Feynman-Kac path integral provides the solution to the diffusion equation. In Sec- 
tion lYLH we report on algorithms and complexity on a classical computer. In the worst case 
for a d- dimensional Feynman-Kac path integral the problem again suffers the curse of dimen- 
sionality which can be broken by Monte Carlo. In Section IVI.2I we indicate the algorithm 
and query complexity on a quantum computer. 

Section PVIIt Eigenvalue approximation. 

One of the most important problems in physics and chemistry is approximating the 
ground state energy governed by a differential equation. Typically, numerical algorithms on 
a classical computer are well known. Our focus is on the Sturm-Liouville eigenvalue (SLE) 
problem where the first complexity results were recently obtained. The SLE equation is also 
called the time-independent Schrodinger equation. In Section lVlI.ll we present an algorithm 
on a classical computer. The worst case query complexity suffers the curse of dimensionality. 
We also state a randomized algorithm. The randomized query complexity is unknown for 
d > 2 and is an important open question. In Section IVII.2I we outline an algorithm for a 
quantum computer. The quantum query complexity is not known when d > 1. It has been 
shown that it is not exponential in d; the problem is tractable on a quantum computer. 

Section IVIII1 Qubit complexity. 

So far we've focused on query complexity. For the foreseeable future the number of 
qubits will be a crucial computational resource. We give a general lower bound on the qubit 
complexity of continuous problems. We show that because of this lower bound there's a 
problem that cannot be solved on a quantum computer but that's easy to solve on a classical 
computer using Monte Carlo. 

A definition of a quantum algorithm is given by flTJ); the queries are deterministic. 
Wozniakowski [77J introduced the quantum setting with randomized queries for continuous 
problems. For path integration there is an exponential reduction in the qubit complexity. 

Section IIX1 Approximation. 

Approximating functions of d variables is a fundamental and generally hard problem. 
The complexity is sensitive to the norm, p, on the class of functions and to several other 
parameters. For example, if p = oo approximation suffers the curse of dimensionality in 
the worst and randomized classical cases. The problem remains intractable in the quantum 
setting. 

Section |Xj Partial Differential Equations. 

Elliptic partial differential equations have many important applications and have been 
extensively studied. In particular, consider an equation of order 2m in d variables. In the 
classical worst case setting the problem is intractable. The classical randomized and quantum 
settings were only recently studied. The conclusion is that the quantum may or may not 
provide a polynomial speedup; it depends on certain parameters. 

Section IXH Ordinary Differential Equations. 

Consider a system of initial value ordinary equations in d variables. Assume that the 
right hand sides satisfy a Holder condition. The problem is tractable even in the worst 
case with the exponent of e~ x depending on the Holder class parameters. The complexity 
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of classical randomized and quantum algorithms have only recently been obtained. The 
quantum setting yields a polynomial speedup. 

Section Kill Gradient estimation. 

Jordan [37] showed that approximating the gradient of a function can be done with a 
single query on a quantum computer although it takes d+ 1 function evaluations on a classical 
computer. A simplified version of Jordan's algorithm is presented. 

Section IXIII1 Simulation of quantum systems on quantum computers. 

There is a large and varied literature on simulation of quantum systems on quantum 
computers. The focus in these papers is typically on the cost of particular classical and 
quantum algorithms without complexity analysis and therefore without speedup results. To 
give the reader a taste of this area we list some sample papers. 

Section IXIV1 Future directions. 

We briefly indicate a number of open questions. 

Ill Overview of Quantum Algorithms 

A quantum algorithm consists of a sequence of unitary transformations applied to an initial 
state. The result of the algorithm is obtained by measuring its final state. The quantum 
model of computation is discussed in detail in [HI [7J El CEO E7J HE] and we summarize it here 
as it applies to continuous problems. 

The initial state |^o) of the algorithm is a unit vector of the Hilbert space 7i u = C 2 (8) 
• • • ® C 2 , v times, for some appropriately chosen integer u, where C 2 is the two dimensional 
space of complex numbers. The dimension of 7i u is 2 V . The number v denotes the number 
of qubits used by the quantum algorithm. 

The final state is also a unit vector of 7i u and is obtained from the initial state \i/jq) 
through a sequence of unitary 2 V x 2 V matrices, i.e., 



The unitary matrix Qf is called a quantum query and is used to provide information to the 
algorithm about a function /. Qf depends on n function evaluations f(ti),...,f(t n ), at 
deterministically chosen points, n < 2 V . The Uq,Ui, . . . ,Ut are unitary matrices that do not 
depend on /. The integer T denotes the number of quantum queries. 

For algorithms solving discrete problems, such as Grover's algorithm for the search of 
an unordered database [2H], the input / is considered to be a Boolean function. The most 
commonly studied quantum query is the bit query. For a Boolean function / : {0, . . . , 2 m — 
1} — > {0, 1}, the bit query is defined by 



Here v = m + 1, \j) G TC m , and \k) G TCi with © denoting addition modulo 2. For a real 
function / the query is constructed by taking the most significant bits of the function / 
evaluated at some points tj. More precisely, as in [27], the bit query for / has the form 



(1) 



Qf\j)\k) = |J>|* ©/(?)>■ 



(2) 



Qf\j)\k) = \j)\k®(3(f(T(j)))), 



(3) 
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where the number of qubits is now v = m' + m" and \j) G 1~L m >, \k) G 7Y m ». The functions 
/3 and r are used to discretize the domain T> and the range 1Z of /, respectively. Therefore, 
(3 : K -f {0, 1, . . . , 2 m " - 1} and r : {0, 1, . . . , 2 m ' - 1} -> V. Hence, we compute / at 
tj = r(j) and then take the m" most significant bits of f(tj) by (3(f(tj)), for the details and 
the possible use of ancillary qubits see [27] . 

At the end of the quantum algorithm, the final state \tpf) is measured. The measurement 
produces one of M outcomes, where M < 2 U . Outcome j G {0, 1, . . . , M — 1} occurs with 
probability Pf(j), which depends on j and the input /. Knowing the outcome j, we classically 
compute the final result of the algorithm. 

In principle, quantum algorithms may have measurements applied between sequences of 
unitary transformations of the form presented above. However, any algorithm with multiple 
measurements can be simulated by a quantum algorithm with only one measurement [8]. 

Let S be a linear or nonlinear operator such that 

S : T -> Q. (4) 

Typically, T is a linear space of real functions of several variables, and Q is a normed linear 
space. We wish to approximate S(f) to within e for / G T. We approximate S(f) using n 
function evaluations • • • , f(t n ) at deterministically and a priori chosen sample points. 

The quantum query Qf encodes this information, and the quantum algorithm obtains this 
information from Qf. 

Without loss of generality, we consider algorithms that approximate S(f) with probability 
p > |. We can boost the success probability of an algorithm to become arbitrarily close to 
one by repeating the algorithm a number of times. The success probability becomes at least 
1 — 5 with a number of repetitions proportional to log<5 -1 . 

The local error of the quantum algorithm ([T|) that computes the approximation <f)f(J), 
for / G T and the outcome j G {0, 1, . . . , M — 1}, is defined by 

e{4> f ,S) = min { a : £ P/(j)>f}, (5) 

L r- \\s(f)-<Pf(j)\\< a } 

where Pf(j) denotes the probability of obtaining outcome j for the function /. The worst 
case error of a quantum algorithm is defined by 

e quant (0,S) = supe^S). (6) 

The query complexity comp query (e, S) of the problem S is the minimal number of queries 
necessary for approximating the solution with accuracy e, i.e., 

comp query (e) = min{T : 3 such that e quant (0, S) < e }. (7) 

Similarly, the qubit complexity of the problem S is the minimal number of qubits necessary 
for approximating the solution with accuracy e, i.e., 

comp qubit (£) = min{ v : 3 such that e quant (0, S) < e }. (8) 
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IV Integration 



Integration is one of the most commonly occurring mathematical problems. One reason is 
that when one seeks the expectation of a continuous process one has to compute an integral. 
Often the integrals are over hundreds or thousands of variables. Path integrals are taken 
over an infinite number of variables. See Section |Vj 

IV. 1 Classical Computer 

We begin with a one dimensional example to illustrate some basic concepts before moving 
to the (i-dimensional case (Number of dimensions and number of variables are used inter- 
changeably.) Our simple example is to approximate 

1(f) = I f{x)dx. 
Jo 

For most integrands we can't use the fundamental theorem of calculus to compute the integral 
analytically; we have to approximate numerically (most real world continuous problems have 
to be approximated numerically). We have to make a promise about /. Assume 

Fi = {/ : [0, 1] -> R | continuous and |/(z)| < 1, x £ [0, 1]} . 

Use queries to compute 

Vf = [f( x l)>---J( x n)\- 

We show that with this promise one cannot guarantee an ^-approximation on a classical 
computer. We use a simple adversary argument. Choose arbitrary numbers Xi,...,x n £ 
[0, 1]. The adversary answers all these queries by answering f(x{) — 0, i = 1, . . . , n. 

What is /? It could be f\ = and f f\(x)dx = or it could be the f 2 shown in 
Figure [H The value of f 2 (x) dx can be arbitrarily close to I. Since fi(x) and fi(x) are 
indistinguishable with these query answers and this promise, it is impossible to guarantee 
an ^-approximation on a classical computer with e < ~. We will return to this example in 
Section Emi 

We move on to the (i-dimensional case. Assume that the integration limits are finite. For 
simplicity we assume we're integrating over the unit cube [0, l] d . So our problem is 

/(/)= f f(x)dx. 

J[0,l] d 

If our promise is that / £ F^, where 

F d = {/ : [0, — > IR. | continuous and \ f{x)\ < 1, x £ [0, l] d ] , 

then we can't compute an ^-approximation regardless of the value of d. Our promise has 
to be stronger. Assume that integrand class has smoothness r. There are various ways to 
define smoothness r for functions of d variables. See, for example, the definition [681 P- 25]. 
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Figure 1: All the function evaluations are equal to zero but the integral is close to one 

(For other definitions see |67j.) With this definition Bakhvalov [I] showed that the query 
complexity is 



This is the worst case query complexity on a classical computer. 

What does this formula imply? If r = the promise is that the functions are only 
continuous but have no smoothness. Then the query complexity is oo, that is, we cannot 
guarantee an ^-approximation. If r and e are fixed the query complexity is an exponential 
function of d. We say the problem is intractable. Following R. Bellman this is also called the 
curse of dimensionality. In particular, let r = 1. Then the promise is that the integrands 
have one derivative and the query complexity is 0(e~ d ). 

The curse of dimensionality is present for many for continuous problems in the worst case 
setting. Breaking the curse is one of the central issues of information-based complexity. For 
high-dimensional integration the curse can be broken for Fd by the Monte Carlo algorithm 
which is a special case of a randomized algorithm. The Monte Carlo algorithm is defined by 



where the Xi are chosen independently from a uniform distribution. Then the expected error 
is 




(9) 




(10) 





where 




is the variance of /. For the promise / 6 the query complexity is given by 
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This is the randomized query complexity on a classical computer. 

Thus Monte Carlo breaks the curse of dimensionality for the integration problem; the 
problem is tractable. Why should picking sample points at random be much better than 
picking them deterministically in the optimal way? This is possible because we've replaced 
the guarantee of the worst case setting by the stochastic assurance of the randomized setting. 
There's no free lunch! 

The reader will note that (ITTj) is a complexity result even though it is the cost of a 
particular algorithm, the Monte Carlo algorithm. The reason is that for integrands satisfying 
this promise Monte Carlo has been proven optimal. It is known that if the integrands are 
smoother Monte Carlo is not optimal; See [68j P- 32]. 

In the worst case setting (deterministic) the query complexity is infinite if / G Fj- In the 
randomized setting the query complexity is independent of d if / G F^. This will provide us 
guidance when we introduce Monte Carlo sampling into the model of quantum computation 
in Section IVIIII 

Generally pseudo-random numbers are used in the implementation of Monte Carlo on 
a classical computer. The quality of a pseudo-random number generator is usually tested 
statistically; see, for example, Knuth [41J. Will (fTTj) still hold if a pseudo-random number 
generator is used? An affirmative answer is given by [70] provided some care is taken in the 
choice of the generator and / is Lipschitz. 

IV. 2 Quantum Computer 

We've seen that Monte Carlo breaks the curse of dimensionality for high- dimensional inte- 
gration on a classical computer and that the query complexity is of order e~ 2 . Can we do 
better on a quantum computer? 

The short answer is yes. Under certain assumptions on the promises, which will be made 
precise below, the quantum query complexity enjoys exponential speedup over the classical 
worst case query complexity and quadratic speedup over the classical randomized query 
complexity. The latter is the same speedup as enjoyed by Grover's search algorithm of an 
unstructured database [26J. 

To show that the quantum query complexity is of order e~ x we have to give matching, 
or close to matching upper and lower bounds. Usually, the upper bound is given by an 
algorithm, the lower bound by a theorem. The upper bound is given by the amplitude 
amplification algorithm of Brassard et al. [12] which we describe briefly. 

The amplitude amplification algorithm of Brassard et al. computes the mean 

1 N ~ l 

SUM(/) = 

i=0 

of a Boolean function / : {0,1,..., N — 1} — > {0,1}, where N = 2 fc , with error e and 
probability at least 8/7T 2 , using a number of queries proportional to min{e _1 , N}. Moreover, 
Nayak and Wu [UJ show that the order of magnitude of the number of queries of this 
algorithm is optimal. Without loss of generality we can assume that e~ l <C N. 
Perhaps the easiest way to understand the algorithm is to consider the operator 

G={i-\i>m)o fi 
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that is used in Grover's search algorithm. Here 

O f \x) = (-l)^ x) \x), xe{0,l} k , 

denotes the query, which is slightly different yet equivalent [H] to the one we defined in (J2]). 
Let 



few 



be the equally weighted superposition of all the states. 

If M denotes the number of assignments for which / has the value 1 then 

M 

SUM(/) = -. 

Without loss of generality 1 < M < N — 1. Consider the space TC spanned by the states 
W = ]l N^M £ W and ^ l) = /J E 

x:/(x)=0 x:/(x) = l 

Then 

|V)=cos(0/2)|Vo)+siii(0/2)|Vi) 
where sin(0/2) = a/ M/N, and 0/2 is the angle between the states |-0) and |V>o)- Thus 



sin 2 I - 



9\ _ M 
2 J ~ iV' 



Now consider the operator G restricted to 7i which has the matrix representation 

q _ f cos $ — sin 6 1 
I sin 9 cos 9 

Its eigenvalues are A± = e ±ld , i = y/—T, and let |£±) denote the corresponding eigenvectors. 
We can express \ip) using the |£±) to get 

|V) = a|£_> + &|£ + >, 

with a, 6 G C, |a| 2 + \b\ 2 = 1. This implies that phase estimation [18] with G p , p = 2\ 
j = 0, ... ,t — 1, and initial state |O)®*| , 0) can be used to approximate either 9 or 2tt — 9 
with error proportional to 2~*. Note that the first t qubits of the initial state determine 
the accuracy of phase estimation. Indeed, let <fi be the result of phase estimation. Since 
sin 2 (#/2) = sin 2 (7r - 9/2), 

v 

0(2- 



V V) M 



with probability at least 8 /it 2 ; see [T21 HE] for the details. Setting t = 6(log£: x ) and 
observing that phase estimation requires a number of applications of G (or queries Of) 
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proportional to e~ x yields the result. (The complexity of quantum algorithms for the average 
case approximation of the Boolean mean has also been studied [33l 152] .) 
For a real function / : {0, 1, . . . , N — 1} — ► [0, 1], we can approximate 

1 N ~ l 

SUM(/) = -X)/(0, 

i=0 

by reducing the problem to the computation of the Boolean mean. One way to derive this 
reduction is to truncate f(i) to the desired number of significant bits, typically, polynomial 
in logs -1 , and then to use the bit representation of the function values to derive a Boolean 
function whose mean is equal to the mean of the truncated real function, see, e.g. [77\. 
The truncation of the function values is formally expressed through the mapping f3 in (J3J). 
Variations of this idea have been used in the literature [3], [271 150] - 

Similarly, one discretizes the domain of a function / : [0, l] d — > [0, 1] using the function r 
in ([3]) and then uses the amplitude amplification algorithm to compute the average 

1 N ^ 

SUM(/) = -^/(x l ), (12) 

i=0 

for Xi E [0, l] d , i = 0,...,N -1. 

The quantum query complexity of computing the average (Tl2|) is of order e _1 . On the 
other hand, the classical deterministic worst case query complexity is proportional to N 
(recall that e^ 1 <C N), and the classical randomized query complexity is proportional to e~ 2 . 

We now turn to the approximation of high-dimensional integrals and outline an algorithm 
for solving this problem. Suppose / : [0, l] d — > [0, 1] is a function for which we are given some 
promise, for instance, that / has smoothness r. The algorithm integrating / with accuracy e 
has two parts. First, using function evaluations f(xi), i = 1, . . . ,n, at deterministic points, 
it approximates / classically, by a function / with error ei, i.e., 

Il/-/ll<ei, 

where || • || is the norm. The complexity of this problem has been extensively studied 
and there are numerous results [HH EJJ, [68] specifying the optimal choice of n and the points 
Xi that must be used to achieve error e\. Thus 

/ fix) dx — I fix) dx + g(x) dx, 
J[o,i] d J[o,i] d J[o,i] d 

where g — f ' — f. Since / is known and depends linearly on the f{xi) the algorithm proceeds 
to integrate it exactly. So it suffices to approximate Jj x , d g(x) dx knowing that \\g\\ < e\. 

The second part of the algorithm approximates the integral of g using the amplitude 
amplification algorithm to compute 

1 N ^ 

SUM(<7) = -5>(^), 

i=0 
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for certain points yi G [0, l] d , with error 82- Once more, there are many results, see [67J 
[68] for surveys specifying the optimal N and the points y^, so that SUM(g) approximates 
Jjo i] d d x with error 62- Finally, the algorithm approximates the original integral by the 
sum of the results of its two parts. The overall error of the algorithm is proportional to 

e = e 1 + e 2 . 

Variations of the integration algorithm we described are known to have optimal query 
complexity for a number of different promises [271 1291 135| [50] . The quantum query com- 
plexity lower bounds for integration are based on the lower bounds of Nayak and Wu [IT] 
for computing the Boolean mean. The quantum algorithms offer an exponential speedup 
over classical deterministic algorithms and a polynomial speedup over classical randomized 
algorithms for the query complexity of high-dimensional integration. The table below sum- 
marizes the query complexity (up to polylog factors) of high-dimensional integration in the 
worst case, randomized and quantum setting for functions belonging to Holder classes F d a 
and Sobolev spaces W^ d . Heinrich obtained most of the quantum query complexity results 
in a series of papers, which we cited earlier. Heinrich summarized his results in [28J where a 
corresponding table showing error bounds can be found. 





Worst case 


Randomized 


Quantum 


w; 4 , 2 < P < 00 
w; id , i< P <2 
wi d 


£ -d/(r+a) 
£ -d/r 
£ -d/r 
£ -d/r 


£ -2d/(2(r+a)+d) 
£ -2d/(2r+d) 
£-pd/ (rp+pd- d) 
£ -d/r 


£ -d/{r+a+d) 
£ -d/(r+d) 
£ -d/(r+d) 
£ -d/(r+d) 



Abrams and Williams [3] were the first to apply the amplitude amplification algorithm 
to high-dimensional integration. Novak [50] was the first to spell out his promises and thus 
obtained the first complexity results for high-dimensional integration. 

V Path Integration 

A path integral is defined as 

1(f) = [ f(x)Kdx), (13) 
Jx 

where /1 is a probability measure on an infinite-dimensional space X. It can be viewed as 
an infinite-dimensional integral. For illustration we give an example due to R. Feynman. In 
classical mechanics a particle at a certain position at time to has a unique trajectory to its 
position at time t\. Quantum mechanically there are an infinite number of possible trajec- 
tories which Feynman called histories, see Figure [2j Feynman summed over the histories. If 
one goes to the limit one gets a path integral. Setting to — 0, t\ — 1 this integral is 

1(f) = [ /(x)Ai(dx), 
Jc[a,i] 

which is a special case of (1131) . 

Path integration occurs in numerous applications including quantum mechanics, quantum 
chemistry, statistical mechanics, and mathematical finance. 
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Figure 2: Different trajectories of a particle 
V.l Classical Computer 

The first complexity analysis of path integration is due to Wasilkowski and Wozniakowski [72J; 
see also [68, Ch. 5]. They studied the deterministic worst case and randomized settings and 
assume that \x is a Gaussian measure; an important special case of a Gaussian measure is a 
Wiener measure. They make the promise that F, the class of integrands, consists of func- 
tions / : X — > K whose r-th Frechet derivatives are continuous and uniformly bounded by 
unity. If r is finite then path integration is intractable. Curbera [20J showed that the worst 
case query complexity is of order e e 13 where (3 is a positive number depending on r. 

Wasilkowski and Wozniakowski [72] also considered the promise that F consists of entire 
functions and that fi is the Wiener measure. Then the query complexity is a polynomial 
in e^ 1 . More precisely, they provide an algorithm for calculating a worst case ^-approximation 
with cost of order e~ p and the problem is tractable with this promise. The exponent p 
depends on the particular Gaussian measure; for the Wiener measure p = 2/3. 

We return to the promise that the smoothness r is finite. Since this problem is intractable 
in the worst case, Wasilkowski and Wozniakowski [72] ask whether settling for a stochastic 
assurance will break the intractability. The obvious approach is to approximate the infinite- 
dimensional integral by a (^-dimensional integral where d may be large (or even huge) since d 
is polynomial in e^ 1 . Then Monte Carlo may be used since its speed of convergence does not 
depend on d. Modulo an assumption that the n-th eigenvalue of the covariance operator of 
fj, does not decrease too fast the randomized query complexity is roughly e~ 2 . Thus Monte 
Carlo is optimal. 

V.2 Quantum Computer 

Just as with finite dimensional integration Monte Carlo makes path integration tractable on a 
classical computer and the query complexity is of order e~ 2 . Again we ask whether we can do 
better on a quantum computer and again the short answer is yes. Under certain assumptions 
on the promises, which will be made precise below, the quantum query complexity is of 
order e~ l . Thus quantum query complexity enjoys exponential speedup over the classical 
worst case query complexity and quadratic speedup over the classical randomized query 
complexity. Again the latter is the same speedup as enjoyed by Grover's search algorithm 
of an unstructured database. 
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The idea for solving path integration on a quantum computer is fairly simple but the 
analysis is not so easy. So we outline the algorithm without considering the details. We start 
with a classical deterministic algorithm that uses an average as in (TT2]) to approximate the 
path integral 1(f) with error e in the worst case. The number of terms N of the this average 
is an exponential function of e . Nevertheless, on a quantum computer we can approximate 
the average, using the amplitude amplification algorithm as we discussed in Section 4.2, with 
cost that depends on log N and is, therefore, a polynomial in e _1 (71] Sec. 6]. 

A summary of the promises and the results in [71] follows. The measure fi is Gaussian 
and the eigenvalues of its covariance operator is of order j~ h , h > 1. For the Wiener measure 
occurring in many applications we have h = 2. The class of integrands consists of functions / 
whose r-th Frechet derivatives are continuous and uniformly bounded by unity. In particular 
assume the integrands are at least Lipschitz. Then 

• Path integration on a quantum computer is tractable. 

• Query complexity on a quantum computer enjoys exponential speedup over the worst 
case and quadratic speedup over the classical randomized query complexity. More 
precisely, the number of quantum queries is at most 4.22e:~ 1 . 

Results on the qubit complexity of path integration will be given in Section IVnH Details 
of an algorithm for computing an ^-approximation to a path integral on a quantum computer 
may be found in [7T1 Sec. 6]. 

VI Feynman-Kac Path Integration 

An important special case of a path integral is a Feynman-Kac path integral. Assume that 
X is the space C of continuous functions and that the measure \i is the Wiener measure w. 
Feynman-Kac path integrals occur in many applications; see |23j. For example consider the 
diffusion equation 



where u G M, t > 0, V is a potential function, and v is an initial condition function. Under 
mild conditions on v and V the solution is given by the Feynman-Kac path integral 



The problem generalizes to the multivariate case by considering the diffusion equation 




z(u, 0) = v(u) 




(14) 




z(u, 0) = v(u), 
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with u G M. d , t > 0, and V,v : M. d — > R, the potential and the initial value function, 
respectively. As usual, A denotes the Laplacian. The solution is given by the Feynman-Kac 
path integral 



where C is the set of continuous functions x : M + — > M. d such that x(0) = 0. 

Note that there are two kinds of dimension here. A Feynman-Kac path integral is infinite 
dimensional since we're integrating over continuous functions. Furthermore u is a function 
of d variables. 

VI. 1 Classical Computer 

We begin with the case when u is a scalar and then move to the case where u is a multivariate 
function. There have been a number of papers on the numerical solution of (jl4jl : see, for 
example [!%] . 

The usual attack is to solve the problem with a stochastic assurance using randomization. 
For simplicity we make the promise that u — 1 and V is four times continuously differentiable. 
Then by Chorin's algorithm [TO], the total cost is of order e~ 2 - b . 

The first complexity analysis may be found in Plaskota et al. |58j where a new algorithm 
is defined which enjoys certain optimality properties. They construct an algorithm which 
computes an ^-approximation at cost of order e~' 25 and show that the worst case complexity 
is of the same order. Hence the exponent of e _1 is an order of magnitude smaller and with 
a worst case rather than a stochastic guarantee. However, the new algorithm requires a 
numerically difficult precomputation which may limit its applicability. 

We next report on multivariate Feynman-Kac path integration. First consider the worst 
case setting with the promise that v and V are r times continuously differentiable with r 
finite. Kwas and Li [33] proved that the query complexity is of order e~ d l r . Therefore in the 
worst case setting the problem suffers the curse of dimensionality. 

The randomized setting was studied by Kwas [32] ■ He showed that the curse of di- 
mensionality is broken by using Monte Carlo using a number of queries of order e~ 2 . The 
number of queries can be further improved to £- 2 /( 1 + 2r / d ) ; which is the optimal number of 
queries, by a bit more complicated algorithm. The randomized algorithms require extensive 
precomputing [321133]. 

VI. 2 Quantum Computer 

Multivariate multivariate Feynman-Kac path integration on a quantum computer was stud- 
ied in [32]- With the promise as in the worst and randomized case, Kwas presents an 
algorithm and complexity analysis. He exhibits a quantum algorithm that uses a number 
of queries of order e~ x that is based on the Monte Carlo algorithm. He shows that the 
query complexity is of order 5~ 1 /( 1+r / d ) and is achieved by a bit more complicated quan- 
tum algorithm. Just as in the randomized case the quantum algorithms require extensive 
precomputing [42l 143]. 
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VII Eigenvalue Approximation 



Eigenvalue problems for differential operators arising in physics and engineering have been 
extensively studied in the literature; see, e.g. [5] HHJ EH [22], [251 BH E3] [65]. Typically the 
mathematical properties of the eigenvalues and the corresponding eigenfunctions are known 
and so are numerical algorithms approximating them on a classical computer. Nevertheless, 
the complexity of approximating eigenvalues in the worst, randomized and quantum settings 
has only recently been addressed for the Sturm-Liouville eigenvalue problem [53] [55] (see 
also pU] for quantum lower bounds with a different kind of query than the one we discuss 
in this article). In some cases we have sharp complexity estimates but there are important 
questions that remain open. 

Most of our discussion here concerns the complexity of approximating the smallest eigen- 
value of a Sturm-Liouville eigenvalue problem. We will conclude this section by briefly 
addressing quantum algorithms for other eigenvalue problems. 

In the physics literature this problem is called the time-independent Schrodinger equation. 
The smallest eigenvalue is the energy of the ground state. In the mathematics literature it 
is called the Sturm-Liouville eigenvalue problem. 

Let Id = [0, l] d and consider the class of functions 



Q = \ q : I d - [0, 1] 



Qq 

q, Djq := — G C(I d ), ||£>jg||oo < 1, ||g||oo < 1 \ , 



dxj 

where || • denotes the supremum norm. For q G Q, define L g := — A + q, where A = 
Ylj=i d 2 /dx'j is the Laplacian, and consider the eigenvalue problem 

h q u = Xu, x G (0,l) d , (15) 
u(x) = 0, x G dl d . (16) 

In the variational form, the smallest eigenvalue X = X(q) of ffT5"l [T6|) is given by 



j Ji E; = i[^)f+# a w dx 

o^ueh^ fj u 2 (x) dx 



A(g)= min J/ ^- u J J1 * , (17) 



where Hq is the space of all functions vanishing on the boundary of Id having square inte- 
grable first order partial derivatives. We consider the complexity of classical and quantum 
algorithms approximating X(q) with error e. 



VII. 1 Classical Computer 

In the worst case we discretize the differential operator on a grid of size h = G(e _1 ) and obtain 
a matrix M e = —A £ + B £ , of size proportional to e~ d x e~ d . The matrix A e is the (2d + 1)- 
point finite difference discretization of the Laplacian. The matrix B £ is a diagonal matrix 
containing evaluations of q at the grid points. The smallest eigenvalue of M £ approximates 
X(q) with error 0(e) [73] [73]. We compute the smallest eigenvalue of M £ using the bisection 
method [22, p. 228]. The resulting algorithm uses a number of queries proportional to s~ d . 
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It turns out that this number of queries is optimal in the worst case, and the problem suffers 
from the curse of dimensionality. 

The query complexity lower bounds are obtained by reducing the eigenvalue problem to 
high-dimensional integration. For this we use the perturbation formula [53, 55f 

X(q) = X(q) + f ( q ( x ) - q( x )) u 2 q (x) dx + O (\\q - q\\l) , (18) 
J Id 

where q, q £ Q and Uq is the normalized eigenfunction corresponding to A(g). 

The same formula is used for lower bounds in the randomized setting. Namely, the query 
complexity is 

Moreover, we can use (TTSI) to derive a randomized algorithm. First we approximate q by 
a function q and then approximate X(q) by approximating the first two terms on the right 
hand side of (ITS]) ; see Papageorgiou and Wozniakowski j55j for d = 1, and Papageorgiou [53] 
for general d. However, this algorithm uses 

Q/ £ -max(2/3,<i/2)\ 

queries. So, it is optimal only when d < 2. Determining the randomized complexity for 
d > 2 and determining if the randomized complexity is an exponential function of d are 
important open questions. 

VII. 2 Quantum Computer 

The perturbation formula ( ITS]) can be used to show that the quantum query complexity is 

As in the randomized case, we can use (ITS]) to derive an algorithm that uses 0(e~ d ^ 2 ) quantum 
queries. The difference between the quantum algorithm and the randomized algorithm is 
that the former uses the amplitude amplification algorithm to approximate the integral on 
the right hand side of ( ITS]) instead of Monte Carlo. The algorithm is optimal only when 
d = 1 plES]. 

For general d the query complexity is not known exactly, we only have the upper bound 
0(e~ p ), p < 6. The best quantum algorithm known is based on phase estimation. In 
particular, we discretize the problem as in the worst case and apply phase estimation to the 
unitary matrix 

where 7 is chosen appropriately so that the resulting phase belongs to [0, 2tt). We use 
a splitting formula to approximate the necessary powers of the matrix exponential. The 
largest eigenvalue of A e is 0(e~ 2 ) and ||qi|oo < 1- This implies that the resulting number of 
queries does not grow exponentially with d. 

Finally, there are a number of papers providing quantum algorithms for eigenvalue ap- 
proximation without carrying out a complete complexity analysis. Abrams and Lloyd [2] 
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have written an influential paper on eigenvalue approximation of a quantum mechanical 
system evolving with a given Hamiltonian. They point out that phase estimation, which 
requires the corresponding eigenvector as part of its initial state, can also be used with an 
approximate eigenvector. Jaksch and Papageorgiou [36J give a quantum algorithm which 
computes a good approximation of the eigenvector at low cost. Their method can be gen- 
erally applied to the solution of continuous Hermitian eigenproblems on a discrete grid. It 
starts with a classically obtained eigenvector for a problem discretized on a coarse grid and 
constructs an approximate eigenvector on a fine grid. 

We describe this algorithm briefly for the case d = 1. Suppose N = 2 k and Nq = 2 k ° are 
the number of points in the fine and coarse grid, respectively. Given the eigenvector for the 
coarse grid \U^ N °^), we approximate the eigenvector for the fine grid \U^) by 

|0} + |l) \ 8Mfl) 
V2 ) 

The effect of this transformation is to replicate the coordinates of \U^N )) 2 k ~ k ° times. The 
resulting approximation is good enough in the sense that the success probability of phase 
estimation with initial state \Uh) is greater than |. 

Szkopek et al. [M] use the algorithm of Jaksch and Papageorgiou in the approximation 
of low order eigenvalues of a differential operator of order 2s in d dimensions. Their paper 
provides an algorithm with cost polynomial in e~ l and generalizes the results of Abrams and 
Lloyd [2]. However, [64] does not carry out a complexity analysis but only considers known 
classical algorithms in the worst case for comparison. 

VIII Qubit Complexity 

For the foreseeable future the number of qubits will be a crucial computational resource. We 
give a general lower bound on the qubit complexity of continuous problems. 
Recall that in (pQ) we defined a quantum algorithm as 

Iff) = U T Q f U T -xQf ■ ■ ■ UxQfUo^o). 

where \tpo) and \ipf) are the initial and final state vectors, respectively. They are column 
vectors of length 2 V . The query Qf, a 2 V x 2 V unitary matrix, depends on the values of / at 
n < 2 V deterministic points. That is 

Qf = Qf(f(x 1 ),...,f(x n )). 

It's important to note that in the standard model the evolution is completely deterministic. 
The only probabilistic element is in the measurement of the final state. 
For the qubit complexity flSJ) we have the following lower bound 

oomp£*(e, S) = fi(logcomp™(25, S)). (19) 

Here S specifies the problem, see 01]), and comp^j blt (£:, S) is the qubit complexity in the 
standard quantum setting. On the right hand side comp^ ry (£:, S) is the query complexity 
on a classical computer in the worst case setting. See [77] for details. 
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We provide an intuition about (|T9|) . Assume 2 V function evaluations are needed to solve 
the problem specified by S to within e. Note that v qubits are needed to generate the Hilbert 
space 7i u of size 2 U to store the evaluations. 

Equation ( JT9l) can be interpreted as a certain limitation of the standard quantum setting, 
which considers queries using function evaluations at deterministic points. We'll show why 
this is a limitation and show we can do better. 

Consider multivariate integration which was discussed in Section HVl We seek to approx- 
imate the solution of 

S(f)= [ f(x)dx. 
Jo,i] d 

Assume our promise is / G F where 

F d = {/ : [0, l] d -> R | continuous and \ f(x)\ < l,x G [0, l} d ] . 

For the moment let d — 1. We showed that with this promise we cannot guarantee an 
e- approximation on a classical computer with e < \. That is, 

By m 

com Pstd blt ( £ ' 5 ) = °°- 

If the qubit complexity is infinite even for d — 1 its certainly infinite for general d. But we 
saw that if classical randomization (Monte Carlo) is used 

Thus we have identified a problem which is easy to solve on a classical computer and is 
impossible to solve on a quantum computer using the standard formulation of a quantum 
algorithm (Q. Our example motivates extending the notion of quantum algorithm by per- 
mitting randomized queries. The quantum setting with randomized queries was introduced 
by Wozniakowski [77]. The idea of using randomized queries is not entirely new. Shor's 
algorithm [60] uses a special kind of randomized query, namely, 

Qx\x) = \jx modiV), 

with j = 0, . . . , N — 1 and a random x from the set {2, 3, . . . , N — 1}. 

In our case, the randomization affects only the selection of sample points and the number 
of queries. It occurs prior to the implementation of the queries and the execution of the 
quantum algorithm. In this extended setting, we define a quantum algorithm as 

= U Tu Qf^UT u -iQ ffi , ■ ■ ■ £/iQ/,u Wo), (20) 
where u is a random variable and 
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and the Xj w are random points. The number of queries and the points chosen 
at random initially and then remain fixed for the remainder of the computation. Note that 
(|20l) is identical to (CD) except for the introduction of randomization. Randomized queries 
require that we modify the criterion ([6]) by which we measure the error of an algorithm. 
One possibility is to consider the expected error and another possibility is to consider the 
probabilistic error with respect to the random variable to. Both cases have been considered 
in the literature [77] but we will avoid the details here because they are rather technical. 

A test of the new setting is whether it buys us anything. We compare the qubit complexity 
of the standard and randomized settings for integration and path integration. 

VIII. 1 Integration 

We make the same promise as above, namely, / 6 F^. 

• Quantum setting with deterministic queries. We remind the reader that 



Therefore, there is infinite improvement in the randomized quantum setting over the 
standard quantum setting. 

VI 1 1. 2 Path integration 

• Quantum setting with deterministic queries. With an appropriate promise it was shown 



Thus, modulo the log factor, the qubit complexity of path integration is a second degree 
polynomial in e^ 1 . That seems pretty good but we probably won't have enough qubits 
for a long time to do new science, especially with error correction. 

• Quantum setting with randomized queries. Then [77] 



Thus there is an exponential improvement in the randomized quantum setting over the 
standard quantum setting. 




• Quantum setting with randomized queries. Then [77] 

com P ™(e) = e^- 1 ) 
comp? a u n bit (e) = e(loge- 1 ). 



P] that 




e^iog^- 1 ). 



comp?r y (£) = 9(e- 1 ) 
comp^ n bit (e) = e(loge- 1 ). 
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As the analogue of (fl9|) we have the following lower bound on the randomized qubit 
complexity [77] 

compos, S) = Q(log comp™ ran (£, S)), (21) 

where comp^ 1 ™ I I ran (e, S) is the query complexity on a classical computer in the randomized 
setting. 

IX Approximation 

Approximating functions of d variables is a fundamental and generally hard problem. Typ- 
ically, the literature considers the approximation of functions that belong to the Sobolev 
space Wp([0, l] d ) in the norm of L q ([0, l] d ). The condition r/d > 1/p ensures that functions 
in WT([0, l] d ) are continuous, which is necessary for function values to be well defined. Thus, 
when p = oo the dimension d can be arbitrarily large while the smoothness r can be fixed, 
which cannot happen when p < oo. 

For p = oo the problem suffers the curse of dimensionality in the worst and the random- 
ized classical cases |HJ [67] . Recently, Heinrich [30] showed that quantum computers do not 
offer any advantage relative to classical computers since the problem remains intractable in 
the quantum setting. 

For different values of the parameters p,q,r,d the classical and quantum complexities are 
also known [30l EHl E7] . In some cases quantum computers can provide a roughly quadratic 
speedup over classical computers, but there are also cases where the classical and quantum 
complexities coincide. The table below summarizes the order of the query complexity (up to 
polylog factors) of approximation in the worst case, randomized and quantum setting, and 
is based on a similar table in [30] describing error bounds. 





Worst case 


Randomized 


Quantum 


1 < p < q < oo, 
r/d > 2/p-2/q 

1 < p < q < oo, 
r/d < 2/p-2/q 

1 < q < p < oo 


£ -dpq/(rpq-d(q~p)) 
£ -dpq/(rpq-d(q-p)) 
£ -d/r 


£ -dpq/(rpq-d(q-p)) 
£ -dpq/(rpq-d(q-p)) 
£ -d/r 


£ -d/r 

£ -dpq/(2rpq-2d(q-p)) 
£ -d/r 



X Elliptic Partial Differential Equations 

Elliptic partial differential equations have many important applications and have been ex- 
tensively studied in the literature, see [75] and the references therein. A simple example is 
the Poisson equation, for which we want to find a function u : Cl — > R, that satisfies 

— Au(x) = f(x), x G Q 
u(x) = 0, x € <9f2, 

where flcK d , 
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More generally we consider elliptic partial differential equations of order 2m on a smooth 
bounded domain Q C M. d with smooth coefficients and homogeneous boundary conditions 
with the right hand side belonging to C r (Q) and the error measured in the L ro norm; see 
[32] for details. 

In the worst case the complexity is proportional to e~ d ^ r [75J and the problem is in- 
tractable. The randomized complexity of this problem was only recently studied along with 
the quantum complexity by Heinrich [3TT [32] . In particular, the randomized query complexity 
(up to polylog factors) is proportional to 

- max{d/(r+2m), 2d/(2r+d)} 
^ j 

and the quantum query complexity is proportional to 

£— maxjd/ (r+2m), d/(r+d)} 

Thus the quantum setting may provide a polynomial speedup over the classical randomized 
setting but not always. Moreover, for fixed m and r and for d > 4m the problem is intractable 
in all three settings. 



XI Ordinary Differential Equations 

In this section we consider the solution of a system of ordinary differential equations with 
initial conditions 

z\t) = f(z(t)), te[a,b], z(a)=r], 

where / : M d — > M. d , z : [a, b] — > W 1 and 77 G M d with f(rf) ^ 0. For the right hand side 
function / = . . . , f d ], where fj : M d — > M, we assume that the fj belong to the Holder 
class Fj ,a , r + a > 1. We seek to compute a bounded function on the interval [a, b] that 
approximates the solution z. 

Kacewicz [3H] studied the classical worst case complexity of this problem and found it to 
be proportional to s^ l ^ r+a \ Recently he also studied the classical randomized and quantum 
complexity of the problem and derived algorithms that yield upper bounds that from the 
lower bounds by only an arbitrarily small positive parameter in the exponent [39]. The 
resulting randomized and quantum complexity bounds (up to polylog factors) are 

Q( e -V(r+a+l/2-7)\ 

and 

^_l/(r+*fl-7)^ 

where 7 G (0, 1) is arbitrarily small, respectively. Observe that the randomized and quantum 
complexities (up to polylog factors) satisfy 

and 
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respectively. Even more recently, Heinrich and Milla [34J showed that the upper bound for 
the randomized complexity holds with 7 = 0, thereby establishing tight upper and lower 
randomized complexity bounds, up to polylog factors. 

Once more, the quantum setting provides a polynomial speedup over the classical setting. 



XII Gradient Estimation 

Approximating the gradient of a function / : M. d —± R with accuracy e requires a minimum 
of d + 1 function evaluations on a classical computer. Jordan [37] shows how this can be 
done using a single query on a quantum computer. 

We present Jordan's algorithm for the special case where the function is a plane passing 
through the origin, i.e., f(xi, . . . ,Xd) = Ylj=i a j x ji an d is uniformly bounded by 1. Then 
V/ = (ai, . . . , dd) T ■ Using a single query and phase kickback we obtain the state 

JV-l N-l 
ViV 31=0 3d=0 

where N is a power of 2. Equivalently, we have 

1 JV-l N-l 

J2---Y,e 2m{ain+ - +adJd) \ji)---\j d ). 



'N d 

iV h=0 j d =o 



This is equal to the state 



AT-l , JV-l 



31=0 v j d =0 

We apply the Fourier transform to each of the d registers and then measure each register in 
the computational basis to obtain mi, . . . , m^. If Oj can be represented with finitely many 
bits and iV is sufficiently large then rrij/N = aj, j = 1, . . . , d. 

For functions with second order partial derivatives not identically equal to zero the anal- 
ysis is more complicated and we refer the reader to [37] for the details. 



XIII Simulation of Quantum Systems on Quantum Com- 
puters 

So far this article has been devoted to work on algorithms and complexity of problems where 
the query and qubit complexities are known or have been studied. In a number of cases, the 
classical complexity of these problems is also known and we know the quantum computing 
speedup. 

The notion that quantum systems could be simulated more efficiently by quantum com- 
puters than by classical computers was first mentioned by Manin |44| . see also [45] . and 
discussed thoroughly by Feynman [2"4"j . 
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There is a large and varied literature on simulation of quantum systems on quantum 
computers. The focus in these papers is typically on the cost of particular quantum and 
classical algorithms without complexity analysis and therefore without speedup results. To 
give the reader a taste of this area we list some sample papers: 

• Berry et al. [9] present an efficient quantum algorithm for simulating the evolution of 
a sparse Hamiltonian. 

• Dawson, Eisert and Osborne |21j introduce a unified formulation of variational methods 
for simulating ground state properties of quantum many-body systems. 

• Morita and Nishimori [16] derive convergence conditions for the quantum annealing 
algorithm. 

• Brown, Clark, Chuang [13] establish limits of quantum simulation when applied to 
specific problems. 

• Chen, Yepez and Cory [15] report on the simulation of Burgers equation as a type-II 
quantum computation. 

• Paredes, Verstraete and Cirac [56] present an algorithm that exploits quantum paral- 
lelism to simulate randomness. 

• Somma et al. [61] discuss what type of physical problems can be efficiently simulated 
on a quantum computer which cannot be simulated on a Turing machine. 

• Yepez [78] presents an efficient algorithm for the many-body three-dimensional Dirac 
equation. 

• Nielsen and Chuang [48J discuss simulation of a variety of quantum systems. 

• Sornborger and Stewart [62J develop higher order methods for simulations. 

• Boghosian and Taylor [H] present algorithms for efficiently simulating quantum me- 
chanical systems. 

• Zalka [79J shows that the time evolution of the wave function of a quantum mechanical 
many particle system can be simulated efficiently 

• Abrams and Lloyd pQ provide fast algorithms for simulating many-body Fermi systems. 

• Wisner [76] provides two quantum many-body problems whose solution is intractable 
on a classical computer. 
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XIV Future Directions 



The reason there is so much interest in quantum computers is to solve important problems 
fast. The belief is that we will be able to solve scientific problems, and in particular quantum 
mechanical systems, which cannot be solved on a classical computer. That is, that quantum 
computation will lead to new science. 

Research consists of two major parts. The first is identification of important scientific 
problems with substantial speedup. The second is the construction of machines with suffi- 
cient number of qubits and long enough decoherence times to solve the problems identified 
in the first part. Abrams and Lloyd [2] have argued that with 50 to 100 qubits we can solve 
interesting classically intractable problems from atomic physics. Of course this does not 
include qubits needed for fault tolerant computation. 

There are numerous important open questions. We will limit ourselves here to some open 
questions regarding the problems discussed in this article. 

1. In Section |V] we reported big wins for the qubit complexity for integration and path 
integration. Are there big wins for other problems? 

2. Are there problems for which we get big wins for query complexity using randomized 
queries? 

3. Are there tradeoffs between the query complexity and the qubit complexity? 

4. What are the classical and quantum complexities of approximating the solution of 
the Schrodinger equation for a many-particle system? How do they it depend on the 
number of particles? What is the quantum speedup? 
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